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Exclusive production of (f> mesons on hydrogen and deuterium targets is studied in 
the HERMES kinematic region 1 < Q 2 < 7 GeV 2 and 3.0 < W < 6.3 GeV. Spin 
density matrix elements and their Q 2 and t dependences are presented. These data 
are consistent with s-channel helicity conservation in exclusive (j> meson production. 
No statistically significant evidence for the contribution of unnatural-parity-exchange 
amplitudes is found. 

1 Introduction 

Electroproduction of neutral vector mesons from nucleons is described by the interaction 
with the nucleon of qq pair created by virtual photon. The reaction e + TV — ► e' + <f) + N' 
in the single-photon approximation is equivalent to 7* + TV — > <f> + TV', in which the 7* 
has negative squared four-momentum transfer Q 2 and polarization parameter e, the ratio 
of fluxes of longitudinal and transverse virtual photons. The spin transfer from the virtual 
photon 7* to the vector meson V is commonly described Q] [2] in terms of Spin Density 
Matrix Elements (SDMEs). The experimentally determined set of 23 SDMEs for the (j> 
meson is presented here in the Schilling- Wolf |T] representation. From measured values of 
the SDMEs one may examine s-Channel Helicity Conservation (SCHC) in the transition 
7* — >V. The fractional contribution of Unnatural-Parity- Exchange (UPE) amplitudes of the 
process 7* + TV — > V + TV in comparison with Natural-Parity-Exchange (NPE) amplitudes 
can also be derived from SDMEs measurements. Natural-parity exchange indicates that the 
interaction is mediated by a particle of "natural" parity (J p = + , 1~, 2 + , ...), while UPE 
amplitudes describe the exchange of a particle of "unnatural" parity (J p = 0~,1 + ,...). 
Since the <f> meson consists mainly of strange quarks with a very small admixture of other 
quark flavors, it is expected that quark-exchange with the nucleon is suppressed and two- (or 
more) gluon exchange dominates. For p° meson production at HERMES kinematics, both 
two-gluon-exchange and quark-exchange contributions of the same order of magnitude [3J H] 
are needed to describe the data. Hence a comparison of SDMEs for <f> and p° mesons allows a 
test of our understanding of these mechanisms. Further information about these final-state 
interactions can be obtained by measuring the phase difference between amplitudes. 

2 Selection of exclusive mesons 

The experiment was performed with longitudinally polarized electron and positron beams 
at an energy of 27.5 GeV using unpolarizcd hydrogen or deuterium gas targets. The decay 
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products of the </> mesons {<f> — ► if+AT - , branching ratio « 49%) were detected in the HER- 
MES spectrometer [5]. The (j> mesons were selected by requiring 0.99 < Mkk < 1-04 
GeV for the invariant mass of the two hadrons (see Fig. [I]). In addition, the identi- 
fication of the decay particles was required. For the data set obtained in 1996-1997, 
the absence of a threshold pion sig- 
nal from the Cherenkov detector 
was required, while for the data 
from 2000-2002 information from 
the RICH [3] detector was used. 
As can be seen from Fig.[T]the non- 
resonant contribution under the <j) 
meson invariant mass peak is neg- 
ligible. Exclusive events were se- 
lected by the requirement: AE = 
(M% - M 2 )/(2M) < 0.6 GeV, 
where M (Mx) is the mass of the 
nucleon (missing mass). Diffrac- 
tive events were selected by the 
constraint —t' < 0.4 GeV 2 . Here, 
t' = t—to, with t being the squared 
four-momentum transfer from vir- 
tual photon to vector meson and 
—to the smallest kinematically al- 
lowed value of — t at fixed Q 2 and 
energy in the 7* A center-of-mass 
system (W). The fraction of semi- 
inclusive deep-inelastic scattering 
events in the chosen region was de- 
termined to be 1.6%. 
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Figure 1: Two-kaon invariant mass distribution fitted 
with a Breit-Wigner function. The fit indicates a neg- 
ligible non-resonant contribution. 



3 Spin Density Matrix Elements 
3.1 Formalism 

The matrix r? ./ is given by r" = tttt Y\ ,< , >' F-. »' , , ., F*, , , where 

AyA v b J X V X V 2M ^X 1 X^X N X N X v X n X 1 X n X^X^ A v ,A n A^A n ' 

^AvA' a 7 a« ^ s tne amplitude of the process 7*(A 7 ) + N(Xn) —* V(Xy) + N (X N ) (symbols 
in parenthesis denote particle hclicities) and is the normalization factor. The amplitude 
F may be decomposed into the sum of a NPE amplitude T XvX * A Xn and a UPE amplitude 
Us \ ' \ \ ■ Since the contribution to the SDMEs from the NPE amplitudes with nucleon 

V N X~{ Xn r 

spin flip are suppressed by a factor of about — i'/(4M 2 ), only the NPE amplitudes with 
X N = Ajv will be discussed here and the shorthand notation T\ v \^ = T Xv i will be used 
below. The nine Hcrmitian matrices , arc defined in PH. The upper index a = 

X-i x y M 

corresponds to the unpolarized transverse photon, a = 1,2 to the two directions of linear 
polarization, a = 3 - circular polarization, a — 4 corresponds to the longitudinal photon and 
a = 5— 8 arc attributable to the interference of longitudinal and transverse photons. Without 
using different lepton beam energies contributions of the transverse and longitudinal photons 
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cannot be disentangled, hence the matrices rfj and rfj cannot be measured separately; only 

their combination r? 4 = r% + erf a is determined. For this case the normalization factor Af 

l j l j L j 

can be found from the condition tr{r 04 } = Tq 4 + rj 4 + r 04 1 _ 1 = 1. Eight matrix elements 
of the matrices rfj, rjj, rfj, refered to as "polarized" SDMEs, measured for the first time 
by HERMES, are presented in Fig. [2] (yellow background). The other 15 matrix elements, 
called "unpolarized," do not require a polarized beam. 

The 3-dimensional angular distribution of the scattered lepton and the decay products 
is described by the following angles (see detailed definitions in [7]): $ is the angle between 
the <p meson production and lepton scattering planes in the 7*p center-of-mass system. The 
polar (p and azimuthal angles of the decay K + are presented in the ^-meson rest frame. 

The extracted SDMEs will be presented based on a hierarchy of NPE helicity amplitudes: 
I Too | ~ |Tu| > |T i| > |Ti | ~ |Ti_i|, established for the first time in Ref. [8] for p° pro- 
duction. This hierarchy is valid at asymptotically high value of Q 2 , and was experimentally 
confirmed for exclusive p° production at HERMES kinematics [4]. The SDMEs are catego- 
rized into five classes according to this hierarchy. The two classes A and B describe only 
SCHC transitions. Classes from C to E contain also spin flip transitions. Class A comprises 
SDMEs with dominant contributions proportional to | Too | 2 or |Tn| 2 - the helicity -conserving 
amplitudes describing the transitions j L — > Vl and 7^ — > Vt- Class B SDMEs correspond to 
the interference of Too and Tn amplitudes. The main terms for the unpolarized (polarized) 
SDMEs are proportional to the real (imaginary) part of TooT^. In fact, as a general rule for 
the classes B - E, the dominant contribution of the unpolarized (polarized) SDMEs is pro- 
portional to the real (imaginary) part of the product of two amplitudes. Class C contains 
SDMEs with dominant terms that are products of the s-channel helicity non-conserving 
amplitude Tu (corresponding to the j T — * Vl transition), and Tq* or T^ (for Tqq the Tn 
contribution is actually quadratic). Classes D and E are composed of SDMEs in which the 
main terms contain a product of the small helicity- flip amplitudes T10 (7^ — > Vt) and T_n 
(7J — ► V-t), respectively, multiplied by T* v 

3.2 Extracted SDMEs 

The SDMEs are obtained by minimizing the difference between the 3-dimensional (cos 0, <f>, $) 
matrix of the data and the analogous matrix from fully reconstructed simulated events. An 
8x8x8 binning was used for the variables cos©, <fi, <&. The simulated events were gener- 
ated with uniform angular distributions and re- weighted in an iterative procedure with the 
angular distribution W(cos 0, (f>, rfA Q], where the spin density matrix elements r" , 

were treated as free parameters. The best fit parameters were determined using a binned 
maximum log-likelihood method. The minimization itself and the error calculation were 
performed by MINUIT. 

The extracted SDMEs for the kinematic region 1 < Q 2 < 7 GeV 2 , 3 < W < 6.3 GcV, 
and < — t' < 0.4 GeV 2 , are presented for p° and <f> meson data in Fig. [2] The average 
kinematic values are (Q 2 ) = 1.9 GeV 2 , (W) = 4.8 GeV and (-*') = 0.13 GcV 2 . The 
experimental uncertainties are larger for the eight polarized SDMEs due to the imperfect 
lepton beam polarization (0.53), and the small kinematic factor \J\ — e ((e) = 0.8) by which 
the polarized SDMEs are multiplied. 

In Fig. [2] the SDMEs are shown multiplied by certain factors to make the coefficients 
of the dominant amplitude products equal to unity. The elements of class A are presented 
in the figure in such a way that their main terms are proportional to |Tn| 2 , in particular 
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HERMES PRELIMINARY 

■ p° proton, <Q 2 >=1.9 GeV 2 , <W>=5 GeV 
• & proton and deuteron 

A: y*l & Yt ~>V° 



B: Interference Yl^V° & Yt -»V° 
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Figure 2: The 23 SDME's extracted for p° production on the proton (squares) and for 
<f> meson production on proton and deuteron (circles). The inner error bars represent the 
statistical uncertainties, while the outer ones indicate the statistical and systematic uncer- 
tainties added in quadrature. The unshaded (shaded) areas indicate unpolarized (polarized) 
SDME's. For easier interpretation the set of SDMEs was divided in five classes (see text). 



1 — Tqo is chosen. The SDMEs of class A are similar for p° and 4>. The elements of class B are 
also close to each other for both vector mesons. For <f> mesons the values of the elements of 
classes C, D and E fluctuate near zero supporting the validity of SCHC. For p° mesons the 
elements of classes C, D, E with significantly non-zero values indicate that there exists also 
a production mechanism which does not conserve s-channel helicity. We note, that small 
violation of SCHC was observed for <j> production by the HI Collaboration [9]. 

The dependences of <fi meson SDMEs on Q 2 and t' are shown in Fig. [3J For the presen- 
tation of the Q 2 (f) dependence, an additional bin of 0.5 < Q 2 < 1 GeV 2 ( 0.4 < — t' < 1 
GeV 2 ) is included. As seen also in Fig. [21 the SDMEs from classes A and B are non-zero 
(in agreement with a GPD model calculation [3]). One can see some dependence on Q 2 for 
these elements, however a t 1 dependence is not observed. We note that the SCHC relations 
r\_ 1 = — lm{r\_{\, Re{rf } = — Im{rf }, Rc{rf } = Im{r£ } are fulfilled in every kine- 
matic bin in Q 2 and t' . The elements from classes C to E, related to spin- flip transitions, 
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HERMES PRELIMINARY 
■ 1996-2000 data 
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Figure 3: The kinematic dependence of <fi SDMEs on Q 2 and t'. Proton and deuteron data 
are combined. The inner error bars represent the statistical uncertainties, while the outer 
ones indicate the statistical and systematic uncertainties added in quadrature. 
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should mainly depend on t'; the present experimental uncertainty is insufficient to observe 
any dependence. 



3.3 Phase Difference of T n and T 00 

The interference between the dominant amplitudes Tn and Too depends on their phase 
difference S which can be evaluated as follows: tan S = [Im{rj' } + Re{rf }]/[R.e{rf } — 
Im{rf }]- This results in 6 = 33.0 ±7.4 deg. for the combined proton and deuteron <j> meson 
data. For the first time, the sign of 5 is determined to be positive. We note that in the 
GPD-based calculations of Ref. [3], the value of 6 is found to be ~ 3 degrees. This appears 
to be inconsistent with the above result. 



3.4 Unnatural-Parity Exchange 

Without assuming SCHC, the hypothesis of natural parity exchange in the t channel alone 
leads to the following relations: U ± = 1- r$ - 2r i-i - 2r ii ~ 2r i-i = 0, Z7 2 = rfi+r^ = 0, 
and U 3 = rfj + rf_ x = 0. 



As presented in Fig. 0] all those re- 
lations are approximately fulfilled, in- 
dicating that UPE amplitudes do not 
contribute significantly to exclusive </> 
meson production. A non-zero signal 
of UPE is evidence for the existence of 
quark-antiquark exchange, which is ob- 
served for exclusive p° (see Fig. @|. 
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Figure 4: Q 2 dependence of the linear combina- 
tions of SDMEs Ui, U 2 , U 3 for cp and p" (only on 
the top panel). Values integrated over the range 
1 < Q 2 < 7 GeV 2 are presented in the first bin 
Q 2 < 1 GeV 2 by open symbols. Only statistical 
errors are presented. 
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